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Abstract 



We have calculated S {q) and the single particle distribution function < n^> 
for N holes in the t — J model on a non-square \/8 x \/32 16-site lattice with 
periodic boundary conditions. We justify the use of this lattice by appealing 
to results obtained from the conventional 4x4 16-site cluster, and an undoped 
32-site system, each having the full square symmetry of the bulk. This new 
cluster has a high density of k points along the diagonal of reciprocal space, 
viz. along k = A;(l, 1). The results clearly demonstrate that when the single 
hole problem has a ground state with a system momentum of = 1), the 
resulting ground state for N holes involves a shift of the peak of the system's 
structure factor away from the antiferromagnetic state q = (tt, tt). This shift 
effectively increases continuously with N. When the single hole problem has 
a ground state with a momentum that is not equal to = (^, ^), something 
that may easily be accomplished through the use of the t — t' — J model with 
t' /t small and positive, then the above-mentioned incommensurability for N 
holes is not found - the maximum of S{q) remains at g = (tt, tt) for all N. The 
results for the incommensurate ground states can be understood in terms of 
rigid-band filling: the effective occupation of the single hole k = (±|,±|) 
states is demonstrated by the evaluation of the single particle momentum dis- 
tribution function < n^>. Unlike previous studies, we show that for the many 
hole ground state the occupied momentum states are indeed k = (±|,ib|) 
states. These conclusions are in agreement with the predictions for the spiral 
phase made by Shraiman and Siggia. Further, our results demonstrate that 
in some instances important results of moderately doped Cu02 planes can be 
predicted from a knowledge of the properties of weakly doped planes. 
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Typeset using REVT^ 



I. INTRODUCTION: 



The Cu02 plane based high-temperature superconductors have anomalous normal state 
properties, and it is probable that a complete theory of the superconducting instability will 
first require a theory of this phase. One part of the normal state puzzle involves the spin 
dynamics, and in La2-xSrxCu04 for x = 0.075, 0.14 and x = 0.15 recent experiments 
have demonstrated the existence of incommensurate magnetic fluctuations. An explanation 
of these results is an outstanding theoretical problem. 

One theoretical starting point for these materials is the strong coupling limit of the 
Hubbard model [§], and it has been argued that the simpler t — J model p adequately 
represents the important low-energy physics of these systems. Then, the question is: do 
the predictions of the normal state properties extracted from the t — J model agree with 
experiment? Here we shall focus on the predictions of the magnetic features of these systems 
that can be made from the strong-coupling limit. Numerous theoretical treatments [P-[TT 



of this problem have indeed suggested that some form of a magnetic instability towards 
an incommensurate phase may arise in this model. Unfortunately, not all of these theories 
agree with one another, so more work is required to clarify the situation. 

As elaborated in a review by Dagotto |]I2[, one avenue by which theorists may scrutinize 
theoretical predictions involves the use of exact diagonalization techniques. This allows 
for the complete determination of all eigenstates of a given system. The limitation of this 
method is simply the rather small systems that can be studied, and thus comparisons of 
theory to experiments on bulk systems can be quite limited. Progress is being made, and 
recent sophisticated techniques have been developed to treat somewhat larger Hilbert spaces. 
One finite-size scaling study |jl3| of a doped t — J model yielded the encouraging result that 
the commonly used 4x4 16-site square lattice has only small finite size effects, at least 
for one hole. 

An exact diagonalization study of the t — J model for a variety of carrier densities was 
performed by Moreo et al. Ijl^. These studies were conducted before the above-mentioned 
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experiments and thus Moreo et al. focussed on a search for an incommensurate phase 

that was stable in the thermodynamic hmit. We now know that the incommensurabihty is 
only found in the spin dynamics, and thus different questions are important. To be specific: 
(i) what kind of incommensurability (if any) is actually found in the t — J model, and (ii) 
what are the underlying states that the carriers occupy when in such a state? 

To make clear the relevance of the latter question, consider that using the t — J model 
Shraiman and Siggia |^ have predicted the development of an incommensurate spiral phase 
as a Cu02 plane is doped away from half filling. Implicit in the development of their 
theory of the spiral phase is the assumption that at very low doping levels the carriers 
approximately exist in momentum states corresponding to the ground state of the one-hole 
problem. (The one-hole problem has been studied using a variety of techniques, and it is 
well known [p!5|-[T7|| that the ground state has a system momentum k = (±f , ±f )• This has 



been confirmed by various exact diagonalization studies |12], including the finite-size scaling 



work |T3| mentioned above.) Thus, the resilience of some form of rigid band filling around 
so-called hole pockets is crucial if the instability suggested by Shraiman and Siggia is to 
be validated. Since recent photoemission work on YBa2Cu30Q ^5 have found partial 



evidence for the hole pocket picture in the low doping regime, this question is clearly very 
important. 

The potential success of Ref. 6 in predicting the magnetic features of the moderately 
doped high superconductors is related to an even bigger and more important question: 
can lessons learned from studying the properties of the weakly doped Cu02 planes, e.g. the 
single hole problem, be used to correctly extrapolate to higher doping levels? We now know 
that at least for small but nonzero doping levels this may be the case: For one hole localized 



by a divalent Sr impurity, theory has predicted the ground state |]I9|-[2T|]. The magnetic 



component of the ground state was identified |^ , and based on the above-mentioned semi- 



classical field theory of Shraiman and Siggia [TO] it was realized that a two-fold degenerate 



non-coplanar spin texture was present. (This spin state may be thought of as that resulting 
from a particular superposition of ferromagnetic bonds in a 2D antiferromagnetic lat- 



tice I^Tl.) Then, experiment showed that such a model [0,^ correctly reproduced the 
zero-temperature antiferromagnetic correlation length. More recently it was demonstrated 
that this is also true for nonzero temperatures |^5|]. Lastly, using La NQR |^ it has recently 



been observed that at very low dopings and low temperatures, coexisting with long-ranged 
antiferromagnetic order is a transverse spin freezing - the temperature at which the freez- 
ing occurs may be analytically predicted using either the semiclassical field theory, or 



accurately predicted numerically using the model employed in Ref. 25. These successes 
suggest that perhaps one might be able to continue to extrapolate the semiclassical theory 
to even higher doping levels, and this possibility, along with the experiments [ij consistent 
with the spiral phase [^], were the initial motivation for this paper. 

Thus, here we will present two correlation functions measured using ground states ob- 
tained from an exact diagonalization study of the t — J model for metallic densities of 
carriers. We focus on the two questions mentioned above, viz. (i) is there any evidence 
that at non-zero doping levels the ground state for the t — J model displays any hints of 
incommensurability, and (ii) if so, which single-particle states are occupied in the incom- 
mensurate phase. Our results will provide some justification for the spiral phase arguments 
made by Shraiman and Siggia, as well as for the similarities of the ground states for the 
weak and moderate doping regimes. Our paper is organized as follows. In §11 we introduce 
the cluster on which the exact diagonalization determination of the ground state was accom- 
plished. We justify the use of this non-standard, non-square lattice by appealing to exact 
diagonalization results obtained on other lattices possessing the full square symmetry of the 
plane. In §111 we describe the magnetic properties of N holes subject to the t — J model 
for this cluster; we focus on the static structure factor, S'(g), and show that as the cluster 
is doped, the peak in S{q) shifts away from the antiferromagnetic wave vector. Then, in 
§IV we consider the t — t' — J model, and demonstrate what happens when the single hole 
ground state has a crystal momentum that is not at = (±f , ±f ): simply, the above men- 
tioned incommensurability is no longer found. In §V we analyse the occupation of the single 
particle momentum states - we show that for the one, two, three, and four hole systems. 



the occupation of the associated momentum states is not unhke the situation that would be 
predicted via rigid band fiUing arguments. Finally, in §VI we discuss our numerical results, 
focusing on comparisons to other theories and previous exact diagonalization studies. 

II. DESCRIPTION AND PROPERTIES OF NON-SQUARE 16-SITE CLUSTER: 

Exact diagonalization is now a familiar technique by which studies of systems with small 
Hilbert spaces can be carried out. For two-dimensional, S = ^ quantum spin systems (in- 
cluding doped quantum antiferromagnets, a strong-coupling model of the high-temperature 
superconductors) the most commonly studied Hilbert space corresponds to a square 4x4 
cluster with periodic boundary conditions. Since a cluster of spins is supposed to represent 
a portion of the bulk of the crystal, it has always been thought to be desirable to ensure that 
the symmetry of the bulk be maintained in the cluster. In this section we will introduce a 
non-square cluster of 16 spins with periodic boundary conditions. Then, we will justify the 
use of this lattice by comparing to results for clusters which have the full square symmetry, 
and, in fact, we will see that some of the unphysical results obtained with the square 16-site 
cluster are eliminated when our non-square cluster is used. 

Figure 1 shows a cluster of sites which represent a small portion of a bulk, two- 
dimensional square lattice; in total, it contains 32 sites. Also, this cluster has the full 
4mm point group symmetry of the bulk lattice (though obviously not the translational pe- 
riodicity of the infinite square lattice). We impose periodic boundary conditions on this 
cluster, and this yields the reciprocal lattice vectors shown in Fig. 2a. 

Our non-square 16-site lattice is also shown in Fig. 1 - it is outlined by the rectangle 
elongated along the (1,1) direction, and may be referred to as a VS x ^/32 lattice. Clearly, 
it has a lower point group symmetry, viz. it only possesses a centre of inversion symmetry. 
Imposing periodic boundary conditions, the reciprocal lattice vectors for this cluster are 

— * 

shown in Fig. 2b - note that due to the lack of square symmetry of this cluster k — {kx, ky) 
is not necessarily equivalent to {ky, k^)- 
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Our motivation for choosing this cluster is two fold. Firstly, we wish to dope this lattice 
and determine whether or not there is any sign of incommensurability in the many-hole 
ground state. If the ordering wave vector shifts continuously (with doping) away from the 
ordering wave vector for the commensurate antiferromagnetic insulator state, viz. q — (tt, tt), 
then we should employ a cluster that has as many reciprocal lattice vectors close to (tt, tt) as 
possible. As seen in Fig. 2b, our non-square 16-site lattice has a multitude of k points along 
the zone diagonal that are close to (tt, tt), viz. k — (tt, tt), (^, ^), (|, |), and k — (|, |). 
Secondly, if any incommensurability is found in our studies, we wish to understand the 
origin of the possible instability that leads to the incommensurate state. Thus, if we are 
going to scrutinize the above mentioned theories, we should not eliminate the proposed 
progenitors of the incommensurability. Here we shall focus on whether or not the holes tend 
to form many-hole wave functions that are essentially constructed from a rigid band filling 
of the associated one-hole ground states. The important (low energy) one hole states are 

— * — * 

k — (tt, 0), and k = (f,f), and as shown in Fig. 2b, our non-square 16-site cluster does 
indeed possess both of these reciprocal lattice vectors. Thus, the 16-site non-square cluster 
shown in Fig. 1 is ideal for our purposes i/its lack of square symmetry does not produce 
any anomalous results; we now show that this is indeed the case. 

A. Behaviour of the Undoped Non— Square 16-Site Cluster: 

We have evaluated the ground state, and first excited state, for the 32 and 16-site clusters 
shown in Fig. 1, as well as for the common 4x4 square cluster, for the antiferromagnetic 
Hamiltonian 

H = J^Si-Sj (2.1) 

<ij> 

when an 5* = I spin is placed at every site of the cluster, and periodic boundary conditions 

— * 

are used. For all three clusters the ground state was a A; = singlet; for the two 16-site 
clusters, the ground state energies per spin were found to be very close to one another: 
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-.7018 for the 4x4 cluster, and -.7085 for the non-standard 16-site cluster. Further, the 

— * 

first excited state for all three clusters was always found to be a /c = (tt, tt) triplet; for the 
two 16-site clusters the mass gap (per site) was found to be very close: .0723 for the 4x4 
cluster, and .0740 for the non-square cluster. 

We are interested in the magnetic structure factor of the doped lattice; thus, we must 
be sure that the non-square 16-site cluster does not yield any anomalous results for this 
quantity. The magnetic structure factor corresponds to 

S{q) = 1 < GSI\ [ Y: e''-^-^ S^i ■ S- ]\GSI > (2.2) 



where \GS'^ > is the ground state for n holes having system momentum k, and N is the total 
number of sites. In Fig. 3 we show this quantity for all three clusters; since the reciprocal 
lattice points do not always overlap, the comparison can only be made at certain points. It 
is clear that all three clusters give the same general features. Further, and most importantly 
to this study, for the q along the zone diagonal, the 32-site cluster and our non-square 16- 
site cluster have very similar static structure factors. For example, for q — (|, |) all three 
clusters have near identical values oi S{q). 

It is apparent from these results that no anomalous features arise when the non-square 
cluster is used for an undoped Heisenberg Hamiltonian; we now consider the doped cluster. 



B. Behaviour of the Doped Non Square 16 Site Cluster: 

We have investigated the t — J model, defined by 

H = -t Y: {4aCja + h.c.)+JY{S,-S,-^-n,n,) (2.3) 

<i,j>a <ij> 

where we choose t = 1 and J = .4, as representative of a Cu02 plane. The operators c|^, Cjo- 
are the creation and annihilation operators, respectively, corresponding to the Hilbert space 
which has been reduced by having had all doubly occupied sites integrated out; the notation 
< i,j > implies that only near-neighbour pairs are summed over. 
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We added a single hole to the half-filled, antiferromagnetic insulator; then, the minimum 
energy state was determined for every allowed system momentum for both sixteen site 
clusters. The results, cast in the form of a "band" structure, are shown in Fig. 4. The 
variation of energy with respect to wave vector is seen to be similar for the two clusters, 
although the band width for the non-square cluster is smaller than for the square cluster. 

One intriguing advantage to the use of the non-square cluster is quickly recognized from 
these results. To be specific, for one hole and the t — J Hamiltonian, use of the 4x4 square 
cluster yields the entirely unphysical result that all states with system momenta k = (vr, 0) 
and (f , f ) (and, of course, those k points related to these by the 4mm square symmetry) are 
degenerate; a proof of this fact may be found elsewhere p8|. This is unfortunate since these 
two states will be non-degenerate in the bulk limit. Further, theory predicts that these two 
states are the two lowest energy states assumed by a single hole. Our non-square 16-site 
cluster is very useful in that it contains all of these k points, and also has a sufficiently 
small Hilbert space such that the one hole states can be accessed, hut there is no artificial 
(geometry-imposed) degeneracy between k = (vr,0) and k = (|, |). From Fig. 4 it is seen 
that the single hole ground state for the non-square cluster is k = (|, |), and (vr, 0) is an 
excited state; this is consistent with the conclusions that have been reached regarding the 
single hole problem ||T5| , p!6[| . (One persistent disadvantage found when using this cluster for 
one hole is that for k^ — ky = ivr, the minimum energy states are degenerate; the same 
phenomenon occurs for the square 16-site cluster. Only the degeneracy along k^ + ky = ±7r 
found in the square 16-site cluster is removed when we use the non-square cluster.) 

Summarizing this section, we have introduced a non-square 16-site cluster with periodic 
boundary conditions. Important reciprocal lattice points are present in this lattice, and in 
comparison to the square 4x4 16-site cluster certain artificial degeneracies are lifted. No 
anomalous results were found for the undoped or singly doped non-square cluster. 
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III. INCOMMENSURABILITY VS. NUMBER OF HOLES: 



We have used exact diagonalization to find the ground state of the t — J Hamiltonian for 



the non-square 16-site cluster for one through four holes |29]; this corresponds to doping 
levels of X = .0625 to x = .25, and covers the experimental range of interest for systems that 
have displayed incommensurability 

The ground states, for t = 1 and J = .4, for = 1 and 2 holes have crystal momenta 
k = ±(f , f ), and ±(7r, ivr), respectively. For 3 holes the ground state is highly degenerate 
at the following reciprocal lattice points: ±(— vr, 0), ±( — ^, |), ±(— |, |), ±(— ^, ^), and 
±(0, vr). For 4 holes the ground state is found to correspond to momenta ±(f , |). 

We have calculated the static structure factor, defined in Eq. (2.2), and our results are 
shown in Fig. 5. The maximum of S{q) occurs at a wave vector which shifts from (tt, tt) (the 
antiferromagnetic wave vector) for one and two holes, to (^, ^) for three holes, to (f , f ) 
for four holes. (Note that for ground states with non-zero crystal momenta, one should 
perform an average over the set of ground state wave functions that are degenerate (due to 
the degeneracy of the ground state with respect to differing k points); here, for 1, 2, and 4 
holes, due to the lack of mirror symmetry about the x and y axes for our non-square cluster, 
this does not change the results that are obtained when performing this average, viz. only k 
and —k are degenerate, and S{q) is insensitive to which of these ground state eigenfunctions 
is used. For 3 holes the same structure factor is obtained for all of the degenerate wave 
vectors.) 

As the second hole is added, all that happens is a reduction of the antiferromagnetic 
correlations - this may also be seen in another correlation function, viz. the relative decrease 
of the near-neighbour spin-spin correlation function < Si ■ Sj >. Then, for three and four 
holes, an essentially continuous shift in peak position occurs; the continuous shift in wave 
vector is seen to mimic the experiments of Cheong et al. For more than four holes, S{q) 
is essentially fiat, indicating the effective loss of magnetic correlations in the heavily doped 
materials. 
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It would be desirable to be able to perform the same search for incommensurability on a 
lattice with a high density of k points around (vr, vr) such that the neighbouring k points were 
along the (1,0) and/or (0,1) direction; this is the direction of the incommensurate shifts 
found experimentally |1[]. However, the only lattice (with a small number of sites, and thus 
appropriate for exact diagonahzation studies of a multiply doped cluster) is that of a ladder 
of width two - this cluster would in no way approximate the bulk lattice, and thus we must 
be content with a search for incommensurabilities along the zone diagonal. Further, it has 
been suggested, in a weak coupling theory, that one cannot reproduce the experimentally 
observed shifts in a one-band model; instead, a three-band model is required to produce the 



necessary nesting [|30|. Even if we had used a three-band model, the nature of our cluster 
still restricts us to the set of k points explored here, and thus we do not believe studies 
of S{q) on finite clusters in the strong coupling limit could yield more information on the 
incommensurability than we have found until the technical obstacles associated with doping 



a 32, or 36-site cluster |T3|J3T| with many holes are overcome - this may never be possible. 
Further, only with such progress could the finite-size scaling be carried out to scrutinize the 
observation that a very weak logarithmic maximum of S{q) exists at the incommensurate 
wave vectors. 

IV. INCOMMENSURABILITY IN THE t - t' - J MODEL: 

We have considered the ground state of the t — t' — J model. This model corresponds to 
the Hamiltonian of Eq. (2.3) augmented with a next nearest neighbour hopping: 

Ht' = -t' E {clc^'a + h.c) (4.1) 
<i,i'>a 

where i' is a next near neighbour to i. The inclusion of this new term has been motivated in 
a variety of ways ; here it is extremely useful in showing how the incommensurability 



demonstrated in the above section is changed when the crystal momentum associated with 
the single hole problem is shifted away from k = ±(f , f ). For t = 1. and J = .4, we have 
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added a small positive f, viz. f = .2, and found the ground state for our non-square 16-site 
cluster. The ground state momentum is no longer at k = ±(|, |), but now is found to be 
aX k = ±(7r, 0). This behaviour is consistent with the band-structure predictions for a hole 
moving in an inert background. 

In Fig. 6 we show S{q) for 1, 2, 3, and 4 holes, for the t — t' — J model, using the model 
parameters given above. It is clearly seen that the maximum of the magnetic structure 
factor is always at the antiferromagnetic wave vector, viz. q = {it, it). As the cluster is 
progressively doped, all that happens is a suppression of the antiferromagnetic correlations 
- no shift of S{q) to wave vectors neighbouring the antiferromagnetic (vr, vr) is found. (As 
mentioned above, this is also found when studying the near-neighbour spin-spin correlation 
function < Si ■ Sj >.) This is in marked contrast to the behaviour found in the above 
section: cf. Fig. 5. This simple demonstration seems to suggest that in the strong-coupling 
limit the formation of an incommensurate phase requires the single hole problem to have its 
ground state momentum equal to A; = (^, |). We now examine the single particle momentum 
distribution function to show why this is so. 



V. MOMENTUM DISTRIBUTION FUNCTIONS: 

In the previous two sections we have displayed results obtained from exact diagonalization 
studies that provide evidence for incommensurate correlations in the strong coupling limit of 
a two-dimensional doped antiferromagnetic insulator when the single hole ground state was 
located at A; = ±(|, |). The question that naturally arises is: why does the one hole state so 
profoundly affect the many hole features? In this section we wish to show that one can also 
use the exact diagonahzation results to suggest the progenitor of this incommensurability, 
and subsequently answer this question via a study of the electron and hole momentum 
distribution functions. 

Our approach is very similar to one employed by Stephan and Horsch , as well as 



that more recently given in a very clear presentation made by Ding - in our work we 
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shall follow the notation of Ding. One defines the electron distribution function by 

< n„{q) >=< > . (5.1) 
Similarly, a hole momentum distribution function can be defined: 

< n^iq) >=< c^J^^ > . (5.2) 

Note that in using this definition, the hole distribution function includes a spin idex, a 
feature induced by the constraint of no double occupancy - this property is explained by 



Ding [^]. We wish to track the electron and hole occupations as our cluster is doped from 
one to four holes. To be specific, we wish to ascertain which electron and hole states are 
occupied as the incommensurability found in the previous section develops. 

If one examines these distribution functions for the t' = ground states discussed above, 
one must overcome more unphysical degeneracies; e.g., the four hole < n„{q) > has de- 
generate values for qx — % = ±7r- Further, the analysis is greatly complicated by the 
degeneracy (with respect to the wave vector) of the many-hole ground states; e.g., the un- 
physical degeneracy of the 3 hole state. This problem for the non-square 16-site cluster 
is unique to the pure t — J model. To remove it one can add the second near-neighbour 
hopping t' introduced in §0 - it is known that this hopping amplitude is of opposite sign 



to that of the near neighbour hopping ||32|. We have chosen t' /t = — .1 for a number of 
reasons: (i) with this addition the degeneracies of the many hole ground states are lifted, 
(ii) the single hole ground state remains at A; = ±(f, f ), and in comparison to the t' = 
system, the ordering of the low energy excited states is not changed, and (iii) the two and 
four hole ground states become = states, a property that one certainly would expect 
a bulk system with an even number of holes to possess. As an example of the usefulness 
of including the second near neighbour hopping, note that when t' = 0, the three hole 
ground state on our non-square 16-site cluster is degenerate at the following wave vectors: 
±(-7r, 0), ±(-X' f ±(-f ' f ±(-i. x)' Then, when t' = -.1 is added, one finds 

that this unphysical degeneracy is lifted and the ground state occurs at ±(f , — f )• We wish 
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to stress that identical conclusions to the ones presented below can be reached for any small 



and negative t' p7|. 



For 1, 2, 3, and 4 holes in the t— t' — J model with t = 1, J = .4, and t' = —.1, on our non- 
square 16-site cluster the ground state is found to occur ai k = ±(f , f ), (0, 0), ±(f , — f ), and 
(0,0), respectively. Then, the electron and momentum distribution functions are as shown 
in Figs. 7, 8, 9, and 10. For an odd number of holes we show both spin components of the 
distribution functions; for an even number we show just one of the two spin components 
for both the electrons and holes (since the up and down spin distributions are equivalent). 
Various sum rules, etc., associated with these numbers are discussed at length by Ding []36[ . 



The positioning of these numbers in the figures corresponds to the allowed wave vectors of 
the non-square 16-site cluster, as shown in Fig. 2b. 

The one hole ground state (which was found in the subspace of the total magnetization 
being ^z) clearly shows the large occupation of the electron states within the antiferromag- 
netic Brillouin zone (defined by + \qy\ = vr) except for electrons and holes at the wave 
vector of the ground state, i.e. k = ±(f , f ); only one of these electron states is found to 



be occupied. This is the same result as was found by Ding |^6[, and shows that one may 
associate the momentum of the ground state, and the momentum of the hole state, to be 
one and the same for one hole even in the strong coupling limit. Also, a comparison of our 
Fig. 7 and Fig. 4 of Ding provides evidence that our non-square 16-site cluster has 
a momentum distribution function that behaves in a similar fashion to that found for the 
square 4x4 cluster. 

The two hole ground state is a A; = state (which was found in the subspace of zero total 
magnetization), and as Fig. 8 shows, the occupied electron states are within the antiferro- 
magnetic Brillouin zone except for states at the wave vector of the one hole ground states, 
i.e. k = ±(|,|). Now, unlike the one hole ground state, only holes occupy these states. 
This is precisely the distribution function that one would expect based on rigid band filling 
arguments: the minimum energy states for one hole are at k = ±(|, |), and now for two 
holes both of these states are occupied by holes. Also, this is a very different conclusion from 
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that reached by Stephan and Horsch ||3§] whose results on a twenty site lattice suggested 
that the single hole problem had little to do with the many hole ground state. This may be 
understood in part because for their lattice the important ±(|,±|) states are not present 
- see the discussion in §|VT[ 

We have used a small negative t' to lift the unphysical one hole degeneracies, and it 
is these degeneracies that caused the difference between our results, shown in Fig. 8, and 
those of Ding ||36|, for two holes; this is a further example of the usefulness of including t' . 



Ding found that the t — J two hole ground state was degenerate aX k = (0, 0), ±(7r, 0), and 
±(0, tt) for the square 4x4 lattice. The inclusion of a small negative t' lifts this degeneracy 
and makes the ground state a = state. Then, an analysis of the electron and hole 
distribution function clearly shows the occupation of the (±|,±|) states, consistent with 
rigid band filling. This is to be compared with the occupation of the ±(7r, 0), ±(0, tt) states 
that Ding found in his A; = 0, t' = ground state. 

The three hole ground state (which was again found in the subspace of the total magne- 
tization being |5) is at A; = ±(f , — f ), and the electron and hole distribution functions are 
shown in Fig. 9. In comparison to the two hole case, we now see that the third hole occupies 
the same momentum state as the crystal momentum of the ground state, while the first two 
holes are still found to occupy the ±(|, |) states. This is again consistent with rigid band 
filling. To display this we have provided the minimum energy states for one hole in the 
t — t' — J model in Table I. Note that the first excited state within the antiferromagnetic 
Brillouin zone is at k = ±(^, — ^), and is thus the state that one would expect the third 
hole to occupy. Figure 9 is a vivid demonstration of the hole pockets that one would expect 
from rigid band filling arguments. 

The four hole ground state (which was found in the subspace of zero magnetization) is 
a = state. It has electron and momentum distribution functions, as displayed in Fig. 
10, very similar in character to those of the fewer hole states. Hole pockets around the four 
momenta ±(|, ±|) are clearly in evidence; some small tendency towards an expansion of the 
pockets to form a closed fermi surface may be seen. This is consistent with the assumption 
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16| that the band structure around the antiferromagnetic zone faces is very flat along the 
zone boundary, but steep towards the k = point. These results are also suggestive of a 
crossover from hole pocket states to a Luttinger liquid although this simple set of data 
from a flnite cluster can in no way be considered to conclusively answer such an important 
question. 

VI. DISCUSSION: 

We have suggested the use of a non-square 16-site cluster which includes all the im- 
portant reciprocal lattice points for the one hole problem, and lifts certain unphysical de- 
generacies. We have shown that no anomalous results are found for this cluster, and have 
doped it with a small number of holes. The magnetic structure factor clearly shows the 
movement of its peak with carrier density reminiscent of an incommensurate phase. Since 
we are only working with a flnite cluster, and are incapable of doing a flnite scaling analysis 
with these results, we cannot be sure whether these correlations survive in the bulk limit, 
but experiments suggest that only dynamical (i.e. short-ranged) correlations remain. We 
will present the dynamic structure factor for this model in a future publication, and this 
will allow for a more direct comparison with experiment. 

We have studied the electron and hole distribution functions for the many hole problem. 
They provide clear evidence of the development of hole pockets near the ground state wave 
vectors of the one hole problem; as the doping increases it seems quite possible that the hole 
pockets disappear, and a Fermi surface obeying Luttinger's theorem results. On the basis 
of an assumption of (i) the strong coupling limit, (ii) the one hole ground state's character, 
viz. that it is a = ±(|,±|) state producing long-ranged dipolar spin distortions, and 
(iii) the existence of such hole pockets, Shraiman and Siggia proposed the presence of an 
incommensurate spiral phase - our results strongly support their theory. 

Our results contrast with earlier studies of these same questions. Firstly, Moreo et al. 



did not flnd robust evidence of incommensurability when the doping level of a square 
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4x4 cluster described by the t — J model. Our Fig. 5 seem to be very direct evidence of 



such an underlying instability. In contrast to our use of S{q), Dagotto et al. has used 
the dynamic structure factor for a variety of hole fillings, and did not find any evidence of 
incommensurability. Clearly, our use of a cluster that includes all of the important reciprocal 
lattice vectors, and a high density of k points around the antiferromagnetic wave vector, has 
allowed us to make a more direct study of this problem. 

The work of Stephan and Horsch has been considered by some [0 to have clearly 



demonstrated that the single hole problem has nothing to do with the higher doping levels 
of interest. To be specific, their two hole work showed a Luttinger Liquid with a clear 
Fermi surface, and no hint of hole pockets. Their work was conducted on a number of 
different clusters. Our work brings into question the absoluteness of these conclusions - we 
have clear evidence of hole pockets, and a knowledge of the single hole ground and excited 
states are found to be all that is necessary to predict the behaviour of the single particle 
momentum distribution functions for many holes. Thus, the question that must be answered 
is: how can two studies using the same technique (exact diagonalization) produce such totally 
different conclusions? We feel that because our cluster has the important k = ±(|,±|) 
momentum states, and it is these states that are required to properly incorporate the dipole- 
dipole interactions associated with the spiral instability of Shraiman and Siggia and, we 
do indeed find an incommensurability in these ground states, consistent with experiment, 
differing clusters lead to different hole-hole interactions, and these interactions must strongly 
depend on the momentum states that the holes occupy. 

To emphasize this latter point, we note that the work of Ding led him to conclude 



that some form of rigid band filling did indeed occur for two holes. He, however, thought 
that the two single particle states that combined to produce the two hole ground state 
were k = ±(7r, 0), ±(0, tt) states - no hole pockets are then produced. We found that when 
the unphysical degeneracy of these reciprocal lattice points and those at the faces of the 
antiferromagnetic Brillouin zone are lifted (using a small negative t'), and thus a different 
form of rigid band filling, one displaying hole pockets, is produced. This is again an example 
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of the strong dependence of the hole-hole interactions on the underlying single hole ground 
and first few excited states, and the subsequent character of the many hole ground states. 

Our results are clearly in support of some form of rigid band filling {e.g., see Fig. 10), 
and thus suggest that knowledge gained from the study of the simpler one hole problem can 
(sometimes) be used to understand instabilities occurring at higher carrier densities; e.g., 
the incommensurate spiral phase P|. This is similar to conclusions reached previously by 



one of us for the very weakly doped insulator [P3 H27[|, and lends credence to studies of other 
aspects of this problem, e.g. transport in the normal state, that were based on an assumption 
of rigid band filling |^T| having begun with a strong-coupling description of doped Cu02 
planes. 
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FIGURES 

FIG. 1. The 32-site cluster; the rectangle outlines the non-square 16-site cluster that we focus 

on in tliis paper; tliis cluster is seen to be half of the 32-site cluster. 

FIG. 2. Reciprocal lattice vectors for (a) the 32-site cluster, and (b) our non-square 16-site 
cluster. 

FIG. 3. A comparison of the magnetic structure factor for the undoped square 4x4 16-site, our 
non-square 16-site, and the square 32-site, clusters. The reciprocal lattice points are as follows: 
r = (0,0), X = (7r,0), and M=(7r,7r). 

FIG. 4. Band structures for one hole in the 4x4 and non-square 16-site clusters; we have used 
t = 1 and J = .4. 

FIG. 5. Magnetic structure factors for one, two, three, and four holes for the t — J model on 
the non-square 16-site cluster. 

FIG. 6. Magnetic structure factors for one, two, three, and four holes for the t — t' — J model 
on the non-square 16-site cluster, with t'/t = .2. 

FIG. 7. Distribution functions for (a) electrons, and (b) holes, for the single hole problem. The 
ground state is degenerate at k = ±(|, and here we show the distribution functions for the 
k = (—§,—§) state. The upper (lower) numbers represent the spin up (down) components. The 
energy parameters of the t — t' — J model are t = 1, = — .1, and J = .4. The square outlines the 
antiferromagnetic Brillouin zone. 

FIG. 8. Distribution functions for electrons and holes, for the two hole problem. The ground 

state is a /c = state. The upper (lower) numbers represent the electrons (holes). 

FIG. 9. Distribution functions for (a) electrons, and (b) holes, for the three hole problem. The 
ground state is degenerate at k = ib(|,— ^), and here we show the distribution functions for the 
k = (f ,— f) state. The upper (lower) numbers represent the spin up (down) components. 
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FIG. 10. Distribution functions for electrons and holes, for the four hole problem. The ground 
state is a = state. The upper (lower) numbers represent the electrons (holes) . 
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